



m 


/M ' 

- 3 


NASA Technical Memorandum 106513 
ICOMP-94- 3; CMOTT-94-1 


' / •»•••:*'•*•** ** ’ 


r L. 


i 

Modeling of Wall-Bounded Complex 

O 

-0 vi 

fH fC* 

f 

m 

CD 

m 


Flows and Free Shear Flows 

in u 

^ C 

m 

o 

o 

m 


Tsan-Hsing Shih and Jiang Zhu 

UJ 

uj x: 
CC U 
UL W 

m 

O CJ 

o 


Institute for Computational Mechanics in Propulsion 

^ wi 
O < Q 



and Center for Modeling of Turbulence and Transition 

CD (A 

z ^ « 



Lewis Research Center 

*-*■ Q *«* 
J J 1 



Cleveland, Ohio 

Ui IL ^ 
Q 

O X 



and 

3E UJ • 

•J < 

Ct* 

^ X < 

?noz 



John L. Lumley 

U W 

Q» 

^0 Q 



Cornell University 

O UJ o 

p-« a a; m 



Ithaca, New York 

i ^ a 

1 D J 

i- a uu ^ 

§ CD u 

< 1 cc o 




C NAS 
W A LL 
SHEA 
Cent 



Prepared for the 

ASME Fluids Engineering Conference, Symposium on Boundary Layer 
and Free Shear Flows 

sponsored by the American Society of Mechanical Engineers 



I 


MODELING OF WALL- BOUND ED COMPLEX FLOWS AND FREE SHEAR FLOWS 


Tsan-Hsing Shih and Jiang Zhu 
Institute for Computational Mechanics in Propulsion 
and Center for Modeling of Turbulence and Transition 
Lewis Research Center 
Cleveland, Ohio 44135 

John L. Lumley 
Cornell University 
Ithaca, New York 14853 


ABSTRACT 

Various wall-bounded flows with complex geometries 
and free shear flows have been studied with a newly 
developed realizable Reynolds stress algebraic equation 
model. The model development is based on the invari- 
ant theory in continuum mechanics. This theory en- 
ables us to formulate a general constitutive relation for 
the Reynolds stresses. Pope (1975) was the first to in- 
troduce this kind of constitutive relation to turbulence 
modeling. In our study, realizability is imposed on the 
truncated constitutive relation to determine the coeffi- 
cients so that, unlike the standard k-s eddy viscosity 
model, the present model will not produce negative nor- 
mal stresses in any situations of rapid distortion. The 
calculations based on the present model have shown en- 
couraging success in modeling complex turbulent flows. 

1. INTRODUCTION 

The present study concentrates on complex turbu- 
lent shear flows which are of great interest in propulsion 
systems. These flows are backward-facing step flows, 
confined coflowing jets, confined swirling coaxial jets, U- 
duct flows find diffuser flows. Most of these flows have 
complex structures. For example, the confined coflow- 
ing jet combines several types of flow structures, such as 
the shear layer, jet, recirculation, separation and reat- 
tachment. Accurate prediction of these flows is of great 
importance for engine design in all its key elements. Tur- 
bulent free shear flows (such as mix i n g layers, planar 
and round jets) have been also studied for the purpose 
of examining the performance of turbulence models in 
different benchmark flows. 

The turbulence model used in this study is a newly 
developed realizable Reynolds stress algebraic equation 
model which is fundamentally different from the tradi- 
tional algebraic Reynolds stress models. The present 
model is developed using the invariance theory in con- 
tinuum mechanics. This theory leads to a general con- 


stitutive relation for the Reynolds stress tensor u^u] in 
terms of the mean deformation rate tensor Uij and the 
turbulent velocity and length scales characterized by the 
turbulent kinetic energy k and its dissipation rate e. 
Pope (1975) applied this kind of constitutive relation to 
Rodi’s algebraic Reynolds stress formulation in conjunc- 
tion with the LRR second order closure model (Laun- 
der et al., 1975) and obtained an explicit algebraic ex- 
pression for the Reynolds stresses for a two-dimensional 
mean flow field. Taulbee (1992) was able to extend this 
method to a general three-dimensional flow. We no- 
tice that in RodPs algebraic Reynolds stress formula- 
tion, some assumed concepts are in general not valid for 
most turbulent shear flows, for example, the assumption 
of constant anisotropy of the Reynolds stresses and ne- 
glect of turbulent transport of second moments. These 
assumptions may bring large errors to turbulence mod- 
eling. In addition, an inappropriate second order closure 
model would also add errors to this type of model. In 
this study, Rodi’s formulation was not used. We directly 
impose realizability on the constitutive relation for the 
Reynolds stresses to determine the coefficients in the 
relation. As a result, a realizable explicit expression 
for the Reynolds stresses is obtained for general three- 
dimensional turbulent flows. Some model constants are 
fine-tuned against a backward-facing step flow and then 
tested in other flows. 

The calculations axe performed with a conservative fi- 
nite volume method (Zhu, 1991). Grid independent and 
low numerical diffusion solutions are obtained by using 
differencing schemes of second-order accuracy on suffi- 
ciently fine grids. For wall-bounded flows, the standard 
wall function approach (Launder and Spalding, 1974) is 
used for wall boundary conditions. The results axe com- 
pared in detail with experimental data for both mean 
and turbulent quantities. Calculations using the stan- 
dard h-e eddy viscosity model are also carried out for 
the purpose of comparison. The comparison shows that 



the present realizable Reynolds stress algebraic equation 
model significantly improves the predictive capability of 
k-s equation based models, expecially for flows involving 
massive separations or strong shear layers. In these sit- 
uations, the standard eddy viscosity model overpredicts 
the eddy viscosity and, hence, fails to accurately pre- 
dict wall shear stress, separation, recirculation, etc. We 
find that the success of the present model in modeling 
the above mentioned complex flows is largely due to its 
effective eddy viscosity formulation which accounts for 
the effect of mean shear rates. According to the present 
model, the effective eddy viscosity will be significantly 
reduced by the mean strain rate and maintained at a 
correct level to mimic the complex flow structures. 

2. TURBULENCE MODEL 

2.1 Constitutive Relation. Constitutive relations 
for the Reynolds stresses were derived by several re- 
searchers (Pope, 1975, Yoshizawa, 1984 and Rubinstein 
and Barton, 1990). Shih and Lumley (1993) used the 
invariant theory in continuum mechanics and the gen- 
eralized Cayley- Hamilton formulations (Rivlin, 1955) to 
derive a more (perhaps the most) general constitutive 
relation for the Reynolds stresses under the assumption 
that the Reynolds stresses are dependent only on the 
mean velocity gradients and the characteristic scales of 
turbulence characterized by the turbulent kinetic energy 
k and its dissipation rate e. This relation is 
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where 


Hi = Ui, k u k! i, n 2 = u itk u itk , n 3 = u iik ul k> 
n 4 = ul k u? k , n 5 = u itk u ltk u^, n 6 = u^u^u^, 
n 7 = u itk u l , k u l 2 tm u? m (2) 


Eq.(l) contains 11 undetermined coefficients which are, 
in general, scalar functions of various invariants of the 
tensors in question, for example, SijSij (strain rate) 
and flijQtj (rotation rate) which are (H 2 4- Hi)/ 2 and 
(II 2 — Hi)/2 respectively. The detailed forms of these 
scalar functions must be determined by other model con- 
straints such as realizability, and by experimental data. 

It is noticed that the standard k-e eddy viscosity 
model corresponds to the first two terms on the right 
hand side of Eq.(l). Both the two-scale direct interac- 
tion approximation approach (Yoshizawa, 1984) and the 
RNG method (Rubinstein and Barton, 1990) also pro- 
vided a relation which is the first five terms on the right 
hand side of Eq.(l). 

In this study, for simplicity we truncate Eq.(l) to its 
quadratic tensorial form which is of the same form as 
those developed by Yoshizawa (1984) and Rubinstein 
and Barton (1990). 


2.2 Realizability. Realizability (Schumann, 1977, 
Lumley, 1978), defined as the requirement of the non- 
negativity of turbulent normal stresses and Schwarz ? in- 
equality between any fluctuating quantities, is a basic 
physical and mathematical principle that the solution 
of any turbulence model equation should obey. It also 
represents a minimal requirement to prevent a turbu- 
lence model from producing unphysical results. In the 
following, this principle will be applied to the truncated 
constitutive relation Eq.(l) to derive constraints on its 
coefficients. 

Let us first consider a two-dimensional pure mean de- 
formation in which the deformation rate tensor contains 
only non-zero diagonal components, i.e., 

Uij =0, if 

In this case, the normal stress uiui can be written as 

U\ui 1 , _ k r _ 1 * 2 
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If we define a time scale ratio of the turbulent to the 
mean strain rate as 77 = S kfe , where 5 = yj2 SijSij, 
the above equation can be written as 

— 3 + a2? 7 + i2^ a4 + + a7 ) 7 ? 2 


Physically, we know that u\Ui will decrease due to the 
stretching by E7i t i. However, by realizability ujui should 
not be driven to negative values. Therefore, we require 
that 
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These physically necessary conditions are called the re- 
alizability conditions. Similar analysis of u 2 u 2 and 113U3 
also leads to the above conditions. In addition, it should 
be mentioned that the above analysis also holds for the 
situation of a three-dimensional pure strain rate. These 
conditions can be satisfied in several ways. Among them 
the simplest way is perhaps the following: 
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Two quantities, the turbulent kinetic energy k and its 
dissipation rate s, remain to be determined in Eq.(3). 
To this end, we use the standard k-e model equations 
which are 
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where £ = Qk/e, £1 = (2Q^O- J ) 1 ^ 2 , £1^ = {Uij - 
Uj t i)/2 -f 4 and represents the rotation of 
the coordinate frame. Ai, A 2 , C r i, C r2 and C T $ will be 
taken as constants and determined by comparing calcu- 
lations with experiments. 

It can be seen from the above analysis that realiz- 
ability cannot be fully satisfied if the model coefficients 
(a 2 -d7) are taken as constant, such as those in the stan- 
dard k-€ model and some anisotropic models, such as 
the model of Speziale (1987). In fact, these models sat- 
isfy realizability only in the weak sense, i.e., they only 
ensure the positivity of the sum of the normal Reynolds 
stresses. For more detailed discussion about model co- 
efficients see Shih et al.(1993). 

2.3 Model Equations. The realizable Reynolds 
stress algebraic equation model can be written as 

2 

UiUj — —kSij — Vt(Ui 3 j + Uj,i) 


C e i = 1.44, C s2 = 1.92, cr K = 1, a s = 1.3 

and the other coefficients are taken as 

C T i = ™4, C T 2 ~ 13, C T 3 = —2, A\ — 5.5, A 2 — 1000. 

These values are calibrated against the backward-facing 
step flow of Driver and Seegmiller (1985) for which a 
complete set of experimental data is available for both 
mean and turbulent quantities and they are also found 
to be appropriate for other complex flows studied in this 
work. 

3. APPLICATIONS 

3.1 Diffuser Flows. Two conical diffuser flows 
were calculated, one with a S° total angle (Trupp et 
aL, 1986) and the other 10° (Fraser, 1958). In both 
cases, the flows undergo strong adverse pressure gradi- 
ents but remain attached. Although the flow configura- 
tion looks simple, it is not easy to calculate this type of 
flow accurately, especially for the boundary layer quan- 
tities. Fig.l shows the variation of calculated and mea- 
sured wall friction coefficient Cf with the axial distance 
xjR 0 (R 0 is the inlet duct radius). It is seen that the 
result of the present model is in good agreement with 
the experimental data, while the standard k-s (SKE) 
model overpredicts Cf along almost the entire length of 
the diffuser. The calculated and measured displacement 
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thickness S* are compared in Fig. 2. The comparison 
shows that the SKE model gives a good prediction in 
the upstream region, but deviates significantly from the 
experiment downstream; the present model prediction is 
good in the whole region. Fig. 3 shows the comparison of 
calculated and measured shape factor H. This is the case 
in which the worst agreement with the measurement has 
been found for both models. Nevertheless, the present 
model still performs considerably better than does the 
SKE model. 

3.2 U-Duct Flow. This case is the experiment 
of Monson et al. (1990) conducted in a 180° planar 
turnaround duct. It features flow with large streamline 
curvature. Calculations are compared to the experiment 
taken at a flow Reynolds number of 10 6 . Fig. 4 shows the 
streamlines computed with the present model. A small 
separation region is found at the bend exit. However, 
the SKE model does not predict the flow separation. 
Fig. 5 shows the comparison of calculated and measured 
Cf along the inner wall. The bend is located between 
21.7<s/H<24.8. Both models are seen to behave in 
the same manner and produce large discrepancies in the 
bend region. The reason for this may partially due to 
the use of the wall function which does not respond to 
the severe pressure gradient. 

3.3 Backward- Facing Step Flows. Two back- 
ward facing step flows, measured by Driver and Seeg- 
miller (1985) and Kim et al. (1978), were calculated. 
The former (DS case) has a smaller and the latter (KKJ 
case) a larger step expansion. The computed and mea- 
sured reattachment points are compared in Table 1. The 
calculated reattachment point from the present model 
agrees well with the experiments. Fig. 6 shows the com- 
parison of the computed and the measured static pres- 
sure coefficient C p along the bottom wall. The SKE 
model is seen to predict a premature pressure rise, which 
is consistent with its underprediction of the reattach- 
ment length, while the present model captures the pres- 
sure rise quite well. Fig. 7 shows the comparisons of pre- 
dicted and measured turbulent stresses uu, vv and uv at 
the location x=2 which is in the recirculation region. In 
the KKJ-case, no reliable experimental data exist for the 
turbulent stresses due to the unsteadiness of the flow. 
However, the experimental data of the DS-case is con- 
sidered more reliable because of the smaller unsteadiness 
of the flow. As compared with the results of the SKE 
model in Fig. 7, it is seen that the anisotropic terms in 


the present model increase uu and decrease vv, leading 
to significant improvements in both uu and vv except in 
the near- wall region. On the other hand, the anisotropic 
terms have little impact on uv. The improvement ob- 
tained by the present model for uv is mainly due to the 
reduction in C ^ by strain rate. 


Table 1. Comparison of the reattachment points 


Case 

measurement 

SKE 

PRESENT 

DS 

6.1 

4.99 

5.82 

KKJ 

7± 0.5 

6.35 

7.35 


3.4 Confined Jets. The general features of con- 
fined jets, the experiments of Barchilon and Curtet 
(1964), are sketched in Fig.8. At the entrance, two 
uniform flows, a jet of larger velocity and an ambient 
stream of smaller velocity, are discharged into a cylin- 
drical duct of diameter D 0 . The inlet flow conditions 
can be characterized by the Craya- Curtet number C t . 
The experiment shows that recirculation occurs when 
Ct <0.96. For a given geometry, recirculation as well as 
adverse pressure gradients can be intensified by reduc- 
ing the value of C t at the entrance. The separation and 
reattachment points of the predicted recirculation bub- 
bles are compared with the experimental data in Fig. 9. 
The experiment indicated that as C t decreases, the sep- 
aration point moves upstream while the reattachment 
point remains practically unchanged. The present model 
captures this feature well and predicts both the separa- 
tion and reattachment points much better than does the 
SKE model. The variation of the pressure coefficient C v 
along the duct wall is shown in Fig.10. The pressure dis- 
tribution is governed by the jet entrainment as well as 
the contraction and expansion of the flow caused by the 
recirculation bubble. The decrease in the ambient veloc- 
ity induced by the entrainment gives rise to an adverse 
pressure gradient, while the contraction of stream l i n es 
produces the opposite effect. These two mechanisms in- 
teract more intensely with each other as C t decreases 
and cause the pressure to vary little in the region up- 
stream of the center of the recirculation bubble. How- 
ever, in the downstream part of the recirculation bubble, 
the deceleration of the flow sets up an adverse pressure 
gradient, the slope of which becomes steeper as Ct de- 
creases. Therefore, the ability to capture the location 
of the recirculation center will have a direct impact on 
the prediction of the pressure. Regarding the compari- 
son between predictions and experiments, it is seen that 
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although both models predict practically the same to- 
tal pressure rises which are in excellent agreement with 
the measurements, the present model captures the steep 
pressure gradients better than does the SKE model for 
all of the Ct values. 

3.5 Confined Swirling Coaxial Jets. This is the 
case experimentally studied by Roback and Johnson 
(1983). Fig. 11 shows the general features of the flow. 
At the inlet, an inner jet and an annular jet are ejected 
into an enlarged duct. Besides an annular recirculation 
bubble due to sudden expansion of the duct, a centerline 
recirculation bubble is created by flow swirling. Fig. 12 
compares the calculation of the centerline velocity with 
the experiment. The negative velocity indicates the cen- 
tral recirculation. It is seen that both models predict the 
strength of central recirculation and the front stagnation 
point quite well, but the present model predicts the rear 
stagnation point much better than does the SKE model. 
Fig. 13 shows the comparison of calculated and measured 
mean velocity profiles at x=5.1cm. Both models give 
reasonably good profiles which are within experimental 
scatter, except for the peak values of the axial and radial 
velocities. Both models have been found to give nearly 
the same results in the downstream region, which can 
also be seen from Fig. 12. 

3.6 Turbulent Free Shear Flows. Calculations 
were also performed for a mixing layer, a plane and 
a round jet. The results shown here are only for the 
jets due to the space limitation. Figs. 14 and 15 show 
the comparisons of the self-similar mean velocity pro- 
files from the model predictions and the various mea- 
surements for the plane and round jets, respectively. In 
Fig. 14, the model predictions are compared with the 
measurement of Gutmark and Wygnanski (1976) for the 
plane jet. The predictions given by both the present 
model and the SKE model agree well with the exper- 
imental data. For the round jet, the comparisons are 
made between the model predictions and the measure- 
ments of Rodi (1975) and are shown in Fig. 15. The 
profile distribution of the mean velocity predicted by 
the present model agrees well with Rodi’s data, while 
the SKE model predicts a faster spreading of the round 
jet into the surroundings and a wider distribution. 

4. CONCLUSION 

A realizable Reynolds stress equation model has been 
applied to calculate both complex wall bounded flows 


and free shear flows. The calculations have been com- 
pared with available experimental data. The compar- 
isons show that the present model provides significant 
improvement over the standard k-e eddy viscosity model 
and that the present model is robust and economical as 
well. This indicates that the present model has good 
potential to be a practical tool in engineering applica- 
tions. 
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Fig. 2 Displacement thickness (Case of Trupp et al.) 



Fig.3 Shape factor (Case of Fraser) 
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Fig. 12 Centerline velocity 


Fig. 15 Seif-similar mean velocity profiles for round jet 
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